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OO 

■ We construct a simple translationally invariant, nearest-neighbor Hamiltonian on a chain of 10- 
' dimensional qudits that makes it possible to realize universal quantum computing without any 

^S) , external control during the computational process. We only require the ability to prepare an initial 

computational basis state which encodes both the quantum circuit and its input. The computa- 
tional process is then carried out by the autonomous Hamiltonian time evolution. After a time 
^ polynomially long in the size of the quantum circuit has passed, the result of the computation is 

, obtained with high probability by measuring a few qudits in the computational basis. 

This result also implies that there cannot exist efficient classical simulation methods for generic 
translationally invariant nearest-neighbor Hamiltonians on qudit chains, unless quantum computers 
can be efficiently simulated by classical computers (or, put in complexity theoretic terms, unless 
BPP=BQP). 

■ I. INTRODUCTION 

> 



One of the most important challenges in quantum information science is to identify quantum systems that can be 
controlled in such a way that they can be used to realize universal quantum computing. The quantum circuit model 
abstracts from the details of concrete physical systems and states that the required elementary control operations are: 
(i) initialization in basis states, (ii) implementation of one and two-qubit gates, and (iii) measurement of single qubits 
QQ ' in basis states. Meanwhile, many other models have been proposed such as measurement-based quantum computing 



QQ , |Tl|, [T^ [l^ , adiabatic quantum computing 0, 0] j or topological quantum computing that reduce or modify 

' the set of elementary control operations. However, the common principle underlying all these models is that the 
computation process is always driven by applying a sequence of control operations. 

Instead, we consider a model that does not require any control during the computational process. This model 
00 ^ consists of a quantum system with a Hamiltonian that makes it possible to realize universal quantum computing by 
. the following protocol: (1) prepare an initial state in the computational basis that encodes both the program and 
^ ■ input, (2) let the Hamiltonian time evolution act undisturbed for a sufficiently long time, and (3) measure a small 
• i-H , subsystem in the computational basis to obtain the result of the computation with high probability. We refer to 
' this model as a Hamiltonian quantum computer and more specifically as a Hamiltonian quantum cellular automaton 
^ , (HQCA) provided that the Hamiltonian acts on qudits that are arranged on some lattice, is invariant with respect 
- - ' to translations along the symmetry axis of the lattice, and contains only finite range interactions. Most natural 
Hamiltonians have these properties, so it is important to construct HQCA that are as close as possible to natural 
interactions. 

Hamiltonian QCA are related to the more usual discrete-time QCA (for further review of the different types of 
quantum cellular automata we refer the reader to [l6j). However, while the evolution of discrete-tim e Q CA proceeds 
in discrete update steps (corresponding to tensor products of local unitary operations, see e.g. the states 

of Hamiltonian QCA change in a continuous way according to the Schrodinger equation (with a time-independent 
Hamiltonian). For this reason, Hamiltonian QCA are also called continuous-time QCA [1^. Also, in the HQCA 
model, all the couplings (interactions) are present all the time, while for the the discrete-time QCA, the execution 
of updates on overlapping cells is synchronized by external control. Therefore, the nearest-neighbor interactions of a 
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HQCA have to include a mechanism that ensures that the logical transformations are carried out in the correct order. 

The motivation to consider Hamiltonian computers is threefold. First, it is a fundamental question in the thermo- 
dynamics of computation how to realize computational processes within a closed physical system. Such Hamiltonian 
computers were presented and discussed by Benioff 4*], Feynman f^, and Margolus Second, Hamiltonian quantum 
cellular automata could lead to new ideas for reducing the set of necessary control operations in current proposals 
for quantum computing by using the inherent computational power of the interactions. HQCA are at one end of the 
spectrum of possible implementations; more realistic perspectives for quantum computing could arise by combining 
this model with more conventional models involving external control operations throughout the computation. Third, 
this model can show the limitations of current and future methods in condensed matter physics for simulating the time 
evolution of translationally invariant systems. If evolving with a certain Hamiltonian can realize universal quantum 
computing, then there cannot exist any classical method for efficiently simulating the corresponding time evolution 
unless classical computers are as powerful as quantum computers (BPP=BQP). 

The first theoretical computational models based on a single time- independent Hamiltonian go back to [1, H, . 
However, these Hamiltonian computers were not explicitly designed for realizing universal quantum computing. Mar- 
golus' model has the attractive feature that it is laid out on a 2-dimensional lattice with translationally invariant, 
finite-range interactions. (In Q it was argued that the part of the Hamiltonian responsible for the synchronization in 
a 1-dimensional variant is close to real interaction in solid states.) However, this scheme does not satisfy the require- 
ment (1) since its initial state has to be prepared in a superposition. Building upon Margolus' idea, a translationally 
invariant Hamiltonian universal for quantum computing even if the initial state is restricted to be a canonical basis 
state was given in This model requires 10- local, finite-range interactions among qubits on a 2-dimensional rect- 
angular lattice wrapped around a cylinder. Subsequently, it was established in [13| that nearest-neighbor interactions 
among qutrits on a 2-dimensional lattice suffice. However, the Hamiltonian of [l3l | is translationally invariant only 
when translated over several lattice sites. A different approach was taken by VoUbrecht and Cirac in [21], showing 
that one can implement universal quantum computation with a translationally invariant, nearest-neighbor Hamilto- 
nian on a chain of 30-dimensional qudits. Also, recently another ID translationally invariant Hamiltonian computer 
construction was given by Kay [25|, using particles with dimension d = 31. 

We present two different simplified HQCA constructions on one-dimensional qudit chains. In both models, we think 
of the qudit chain as composed of two registers, data and program. The work qubits we compute on are located at 
a static location in the data register. Driven by the autonomous Hamiltonian time evolution, the program sequence 
contained in the program register moves past the work qubits and the gates are applied to them. After we let the 
system evolve for a time not larger than a polynomial in the length of the program, we measure one or two qudits in 
the computational basis to read out the output of the computation with high probability. 

Our first construction is for a chain of 10-dimensional qudits and is related to the ideas of [2l|. The mechanism 
behind the progress of the program sequence in this particular model can be thought of as the diffusion of a system 
of free fermions on a line. Concurrently with our preprint. Chase and Landahl [23| found another ID construction 
with particles with dimension d = 8. However, their Hamiltonian is not translationally invariant. If we release the 
translational invariance requirement in our d = 10 model, we obtain a d = 8 construction as well. Whether the 
required dimension d = 10 can be decreased while keeping translational invariance in our HQCA model remains an 
open question. 

Our second construction uses qudits with dimension d = 20 and is inspired by [l3|, utilizing a technique of ^ to 
transport the program. Here, the mechanism for the progress of the computation can be thought of as a quantum 
walk on a line. 

The paper is organized as follows. First, in Section fll Al we present the HQCA construction with cell size d = 10 
and analyze the required run-time of this model in Section IIIBI Second, we give the HQCA construction with cell 
size d = 20 in Section [HI Al and discuss the readout procedure and the required run-time in Section [ill Bl We provide 
some useful results for the continuous time quantum walk on a line in Appendix]^ and prove a lemma concerning the 
diffusion of free fermions on a line in Appendix [Bl 



II. THE d = 10 HAMILTONIAN QUANTUM CELLULAR AUTOMATON 

We present a simple universal HQCA on a chain of qudits with dimension d — 10. First, we encode the progression 
of a quantum circuit U on N qubits into a set of states \ip^) of a chain of qudits with length L — poly{N). Second, 
we give a translationally invariant nearest-neighbor Hamiltonian on this chain of qudits, which induces a quantum 
walk on the set of states \^Pa)- Finally, using a mapping to a system of free fermions in ID, we prove that when 
we initialize the qudit chain in an easily determined computational basis state and let the system evolve for a time 
T < Tio = 0{L log L) chosen uniformly at random, we can read out the result of the quantum circuit U with probability 



FIG. 1: a) A quantum circuit consisting of two rounds of gates acting on nearest neighbors, b) The previous circuit with a 
third round of identity gates added. 



Pio ^ I ~ O (^i^^lx) measuring one of the qudits in the computational basis. We then show that this is enough to 
ensure universaHty of our HQCA for the class BQP. 



The Construction 



Encoding a quantum circuit 



The gate set {Toffoli, Hadamard} is universal for quantum computation [2Q|. With only polynomial overhead, one 
can simulate a circuit consisting of these gates using only the gate W (controlled ^ rotation about the y-axis) 
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if it can be applied to any pair of qubits. Let us consider implementing universal quantum computation on a qubit 
chain using only nearest neighbor gates. Let us also restrict the use of the W gate so that the control qubit has to be to 
the left of the target qubit. Using only polynomially many additional swap gates S, one can still do universal quantum 
computation on a qubit chain. Thus given a quantum circuit [/' on N' qubits with poly{N') generic two-qubit gates, 
we can transform it into a circuit [/ on a chain of A'^ = poly{N') qubits with nearest neighbor gates W (with control on 
the left) and S without loss of universality. We then add identity to our gate set and further transform the circuit U 
to have the following form (see Figure[T]). Rewrite the circuit as K rounds of nearest neighbor gates Uk,g G {W, S, /}, 



where gate U^^g belongs to the A:-th round and acts on the pair of qubits Wg, w 



g, UJg+1- 



U ^{UK,N-l.■■UK^)■■■{Ul,N- 



■Ul,l). 



(2) 



We wish to encode the progression of the circuit U into the states of a chain of qudits with dimension 10, with length 
L = poly{N). The basis states of each qudit \q) — \p) (g) \d) are constructed as a tensor product of a 5-dimensional 
program register and a 2-dimensional data register, where p €{.,►, W, 5, /} and d € {0, 1} label the corresponding 
basis states. We start by writing the initial product state 



with pj and dj as in Figure [2] (here we give an example for the circuit in Figure [TJi) : 
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(4) 



The qudit chain has length 



L = 2M = 2KN, 



(5) 



where K is the number of the rounds of gates in ([2]) . The left half of the top (program) register contains K pointer 
symbols ► at positions kN for k = 1 . . . K and empty symbols . everywhere else. The right half holds the program 
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FIG. 2: The initial state \ip) of the qudit chain for the circuit in Figure [T^). 



in the form 
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2,1 ■ 



2.N-1 



I.-.IUk- 



Uk,n-i 



(6) 



first round of gates second round of gates 



last round of gates 



with Uk.g € {W^ S*, /} and each round of gates is preceded by an identity gate. Later we will show a way to execute 
the program while moving it to the left above the data. The bottom (data) register contains N work qubits labeled 
Wn in a canonical basis product state at positions M + n for n — 1 ... TV. We designate as the readout qubit. 
Finally, the rest of the data register (positions 1, . . . , M and M + N + 1, . . . , 2AI) is filled with extra qubits in the 
state |0). 

We now describe how to generate the set of states {|</5cr)} encoding the progression of the quantum computation U. 
Starting from the initial state \ip) we obtain each state \ipa) by applying a particular sequence a of the following 
two update rules. Each of these rules can be applied to two neighboring qudits in the chain. The first rule is: if there 
is an empty spot . to the left of a symbol A e {W, S, 1} (from now on we call it a 'gate') somewhere in the program 
register, move A one step to the left, leaving the data register unchanged. 
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The second rule says: at a place in the qudit chain where a gate A e 
positions and apply the gate A to the qubits in the data register below. 
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(7) 



{W,S,I} meets a 'pointer' switch their 



(8) 



There is only one place where one of these rules can be applied to the initial state ([4]) - the second rule can be applied 
at the place where ► stands to the left of /. After this first step, there are several ways to proceed, generating many 
different states \(pa)- 

Let us look more closely at what happens to the initial state \ip) as we apply some sequence of rules. There are 
K pointers ► in one for each round of gates in the circuit (see Figure [2]). Observe that \ip) is constructed in 
such a way that as a gate Uk,g from the k-th round of gates moves to the left, it meets the k-th (counting from the 
right) pointer ► exactly above the work qubits Wg,Wg+i (see Figure[3fc for a depiction of the second pointer meeting 
the second round of gates). As noted before this is the pair of qubits the gate Uk,g should act on. It could then 
seem that after many applications of the update rules the program moves to the left of the work qubits while the 
corresponding states we obtain have the computation executed on their work qubits as planned. However, before 
drawing that conclusion, we need to consider what happens when a gate Uk^g from the k-th round meets a pointer 
while not above the work qubits (see Figure [3]). First, if this happens above a pair of (extra) data qubits in the 
state 1 00), the qubits stay unchanged after we apply (H]), as we have Uk,g |00) = |00) for the three possible Uk,g (the 
controlled gate W, a swap or the identity). The second possibility is that a gate meets a pointer above the boundary 
of the work qubits, i.e. either above dM,wi or above wn ,dM+N+i- Observe that the identity gates we inserted in 
front of each round of gates in © are the only ones for which this can happen. The extra qubit and the work qubit 
involved then stay unchanged after the application of The extra qubits in the data register thus always remain 
in the state |0) and the only processing in the data register happens on the work qubits. The gates in ^ are applied 
sequentially from left to right. Therefore, when all the gates move to the left half of the chain, the corresponding 
state |(^o-) then contains the result of the quantum circuit U in the state of the work qubits. 



5 



b) 



► II' 



® ® © ® ® © © © © (0) ® (0) 



® ® ® ® ® ® ©©©® ® 



®®®®®®©©©®®® 

FIG. 3: Analysis of the second rule ((Sjl. a) A gate meeting a pointer above two extra qubits in the state |00). b) The gate 
meeting a pointer above the right boundary of the work qubits is the identity gate, c) The only gate meeting a pointer above 
the left boundary of the work qubits is the identity gate. 



2. The Hamiltonian 



We now construct a Hamiltonian whose transition rules will be ([71 13) and their inverses 
It 
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(9) 



Note that the inverse of the rule ^ involves applying A^ to the data qubits, uncomputing the gate. Our Hamiltonian 
is a sum of translationally invariant terms 



L-l 



^io = -E(^+^')(,..+i)' 

where R corresponds to the rules ([7l [8]) and acts on two neighboring qudits as 



i?= > \\A 

Ae{w,s,i} 



(10) 



(11) 



where p stands for the program register and d for the data register of the respective qudit. Similarly, corresponds 
to the inverse rules given in 



3. The computational subspace 



To show that we can obtain a result of a quantum computation by time evolving the state ^ with this Hamiltonian, 
we return to the set of states {|<^ct)} generated using the rules ([71 [5]). There are only finitely many possible sequences 
a. However, when we expand the ruleset by adding the inverse rules the number of possible sequences a of 
rule (forward and backward) applications then becomes infinite. For example, we get the possibility of returning 
back to the initial state \(p), completely undoing the computation. Nevertheless, the space of states {If^So-)} remains 
finite-dimensional for the following three reasons. 

First, the rules ([71 [H [9]) do not change the order of the gates in the program register. Second, once we know the 
locations of the gates, the rest of the program register is uniquely determined. It contains the sequence 




N-l 



repeatedly written out K times from left to right at the positions unoccupied by the gates. Third, we constructed 
the initial state and the rules so that the gates in ([6]) are applied (and uncomputed) to the work qubits sequentially, 
while all the extra qubits remain in the state |0). Given the positions of the gates in the program register of 
we can simply determine how many of the gates Uk^g in ^ have presently been applied to the work qubits. For any 
|</3o-), there can be at most one pointer ► directly above the work qubits wi . . .w^. 
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1 . If there is no pointer located above the work qubits (see Figure [H [3^ and [3j;) , label kp the number of pointers 
already to the right of w^r. The first Nkp gates then have already been applied. 

2. If there is a pointer above the work qubits (see Figure [3Jd), all the gates to the left of this pointer have already 
been applied. 

Let us look at a few examples of the first case. In Figure [21 fcp = and thus no gates have been applied yet. In 
Figures [3^) and[3t) hp = 1 and because iV = 3, we know that the first 3 gates {IWS) have been executed. The state 
of the work qubits is thus equal to 5*2314^12 \w1W2W3), in both Figure [3^) and[5t). 

Therefore, regardless of the particular sequence of update rules (and their inverses) through which we obtained 
\ipa), the state of the data register is uniquely determined by the positions of the gates in the program register. In 
fact, all the states \(pa) with the same position of gates in the program register are the same state. Let us then label 
this state \(pc), where C is a weight M bit string of length L. The M ones in the string C are located at positions 

{aP,...,ai?)}, (13) 

corresponding to the positions of the AI gates in the program register of the state \(pa-)- The number of different 
states in the set {|vc)} is (j(^). Note also that the states \(pc) with different C are mutually orthogonal, as their 
corresponding program registers are in mutually orthogonal product states. The set {If/Jc)} thus forms a basis of an 
(^) -dimensional subspace of the Hilbert space of the qudit chain. 

The transition rules in the Hamiltonian (fTO)) correspond to (l7l[8l[9]). The time evolution of the initial state \(p) (|4]) 
will thus necessarily happen within 

Ti-cornp ^ span{\ipc)}, (14) 

the subspace spanned by the states \ fc)- We call Hcomp the computational subspace. In Section HlB 1\ we show 
that our Hamiltonian pop restricted to the computational subspace has a particularly simple form corresponding to 
a Hamiltonian of a fermionic system. This will allow us to analyze the time evolution of the initial state. 



B. Required Time Analysis 

Our model of computation with the Hamiltonian (fTO)) consists of initializing the qudit chain in the state \(p) (|4]) 
and evolving the system for a time r randomly chosen between and tiq. In this Section, we will show that when we 
measure the output qubit wpf in the data register at time r, we will read out the result of the quantum computation 
U with high probability. Moreover, we will show that the upper bound on the running time of our model, tiq, is 
polynomial in L. 

Each state \ipc) described above can be written as 

IV'c) — \'^c)program ® \Oc)work ® \a)extra, (15) 

where \Trc)program labcls a state of the program register of the chain, the corresponding state of the work qubits is 
\0c)work and the extra data qubits are in the state — \0)f}fra~^^ ■ Note that while each program register 

state \Trc)program is a product state, the corresponding state of the work qubits \Oc)work is entangled. 
The time evolved state \^{t)) is a superposition of the states |<pc)- 

|<^(r)) = e"'-^!"^ \ip) = ^ Cc{t) |<^c) = ^ Cc{t) \nc)program <^ \dc)work <^ \a) extra-, (16) 
C C 

where the sum over C is over all weight M length L bit strings and cc{t) is the amplitude of the state \ fc) in \tf{T)). 

When all of the gates in \'^c) program have moved to the left of the work qubits, the state of the work qubits holds 
the output state of the computation U . For our original setup ([3]), there is only one such state: 

IV'IIIIIIOOOOOO) = / . . ► . . >) program ®U\wi...WN)^,rk®\^) extra- (1^) 

For our model of computation, this would mean that we would need a single amplitude ci...io...o(''') to be large at a 
random time r, which is unlikely. However, a simple modification will turn the tables. Let us pad the qudit chain by 
(/ — 1)M extra sites on the left and M extra sites on the right, with / a constant to be determined later. Let the 
new initial state be 



®(f~l)M 



(18) 
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where is the original initial state ([4]). The original chain had length 2M, so the length of this padded chain is 
L = (/ + 2)M, while the number of gates becomes 2M . As in Section Hi A 3[ let us now label {Iv'n)} the set of states 
which we can obtain from (fT8|) by applying rules ([7l[8l[9|). This time, the label 13 is a weight 2M length L = (/ + 2)M 
bit string, and {aj'°\ . . . , a^^j} denote the positions of the 2M gates in the program register of \^d)- 

The 2M gates in the program register of the initial state \LpP"''^) — |(/3o...oi...i) are located at positions am ' °^ '^^ = 
fM + in with TO = 1 . . . 2M. Only the first M gates are relevant for the computation and the other M are the identity 
gates we inserted as padding. These identity gates leave the data register intact, regardless of whether they move to 
the left of the work qubits or stay to the right of them. Therefore, every state \(Pd) in which the first M (relevant) 
gates have already moved to the left of the work qubits contains the finished computation in the state of its work 
qubits. Let us label the corresponding set of weight 2Af length L = (/ + 2)M bit strings 

Vdone = {D : aif < /A/}. (19) 

For all states |<y5D) with D G 'Ddone^ the state of the work qubits is the same and equal to \9u) — U \wi . . . wn)- Let 
us write the time evolved state of the modified chain using the notation (|15p as 

C_d(t)|v3£)) + ^ CD{T)\TrD)program\ <^\6u)work<^\a)extra (20) 

= Vl -PlO Wfail{T)) + VPIO kdone(T)) ® \Ou)work ® \a)extra- (21) 

The state of the work qubits with the computation done \Qu)work factorizes out in the second term. The two terms 
are orthogonal, as the states of the program register of \ fD) with different D are mutually orthogonal. We will now 
analyze the amplitude of the second term. Assume we have a measurement discriminating different states \t^d) of the 
program register. In that case pio denotes the probability of obtaining \i:d) with D G 'Ddone when measuring at a 
random time r < tiq. 

1. Mapping onto a fermionic system 

We will now prove that when we choose the time r uniformly at random in (O,rio), with tiq = poly{M), the 

probability of finding a state with the computation executed (with om < fM) ispio > 1"*^ (•^) ~ (/ + 2)Af. 

For this, we will analyze the time evolution of ([T5]) with Hiq (fTU]) . We will do it by mapping the states of the 

qudit chain onto states of a chain of spin-i particles, where the positions of gates in the qudit chain will correspond 
to fermionic excitations in the spin chain. By analyzing the latter system, we will bound piQ. 

As in III A3[ the time evolution of the padded initial state (US]) happens only in the computational subspace 

nr,^^ = span{\^n)}. (22) 
The basis states of the computational subspace have the form 



\(PD) = \7TD, 



program 



® ([/(^) . . . WN)^,rk) ® l«)e.t™ : (23) 



where the label D is a weight 2M bit string of length L = {f + 2)M . The states \hd) of the program register are 
canonical basis product states with 

/C/14 . . . Ui.N-i m.i . . . f/2,7v-i • • • IUk,i ■ ■ ■ Uk.n-i (24) 

M 

written out from left to right at the positions {a^^\ . . . , a2M} '^^ ones in D, and the rest of the program register 
is filled with 

^ • ► (25) 



N-l 



written out fK times from left to right. Finally, C/*^^^ stands for the part of the sequence of unitaries (|24p that has 
already been applied to the work qubits when the gates are in positions {ai^\ . . . , a2M}- described the procedure 
to determine J/^^-' from D in Section fll A 31 
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We will now focus on the restriction of Hiq to Ticomp' time-evolved state ([18]) never leaves the computational 

subspace. This restriction has a simple form — it is the negative of an adjacency matrix, where two states are connected 
if one can be obtained from the other by one of the rules ([7l[8l[9|). More exphcitly, two of the states \ ipd) are connected, 
if one can obtain the other by having one gate hop to the left (or right) and applying the appropriate unitary (or do 
nothing) to the work qubits. Consider now 

L-l 

= -J2(m (Oil + |01) (10|),-,+i, (26) 

the Hamiltonian of a line of L = {f + 2)M spin-i particles with a simple hopping interaction, restricted to the 
subspace spanned by states with 2M spin up particles. Observe that with the mapping 

l'/'-D)qudit chain ^ I ^) spin chain ' (^7) 

where D is a weight 2M bit string of length L, the restriction of iJio (|10p to H^^^p (|22p is matrix-wise the same as 
the restriction of Hq to the subspace with 2M up spins, i.e. 

Hiol^pad^^ = i?(}|2Af up spins- (28) 

To show that piQ in (j21[) is large, it will be convenient to analyze the time evolution of the initial state of the spin 
chain 

(29) 

with the spin Hamiltonian (pS)) instead of the time evolution of with JTU]). First, we use the usual Wigner- Jordan 
transformation to define the operators 

b] = aJ...a|_i®|l)(0|^.®I,+i,...,L, (30) 
b, = aJ...a|_i®|0)(l|^.®I,+i,...,L. (31) 

We invite the reader to verify the properties {6^,6^} = 6ijl and — b^^^ — for foj and bj. Thus, foj and bj can be 
viewed as the creation and annihilation operators for a fermion at site j. Rewriting (|26p in terms of (|30p and (|3ip . 
we obtain 

L-l 

Hf^-Y, b]b,+i + h.c, (32) 

a Hamiltonian for a system of free fermions in second quantization. Following our mapping, the initial state [tp^""*) of 
the qudit chain thus corresponds to the state of the fermionic system l^*) = b'^jj^i^i ■ ■ ■ &(y+2)M 1^) with 2M fermions 
on the right end of the line (here |0) is the state with no fermions). We now use the following Lemma (proved in 
Appendix iB)) : 

Lemma 3. Consider the state \^) of 2M fermions on the right end of a line with L — [f + 2)M sites. Let the 
system evolve for a time chosen uniformly at random between and tio with the Hamiltonian given by (j32p and 
measure the number of fermions in the region 1 < x < fM. The probability to measure a number greater than M is 

Pio>^-0(4). 

Let us choose / = 22 and rio — 0{L\ogL) — 0{M log M). Following the mapping we did from our qudit chain 
backwards, this implies that when we initialize the qudit chain of length L = 2AM in \LpP°-'^) as in (fT8|) and let it evolve 
with HiQ for a random time t < tiq , the probability for the chain to be in a state where the gate particles have 
moved sufficiently to the left for the computation to be done (qm < fM) is 



Using this bound on the probability of having the computation done we now proceed to prove universality of our 
model for the class BQP. 
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2. Universality for BQP 



To prove the universality of our model for recognizing languages in the class BQP, let us recall its definition. 
Consider a language L in BQP, a uniform family of circuits U and a problem instance x. When x £ Ly^g, the set of 
instances with the answer 'yes', the probability of the circuit U outputting 'yes' is not smaller than |. On the other 
hand, when x ^ i^es, the probability of the circuit outputting 'yes' is not greater than i. Let us assume the worst 
case for our circuit C/, i.e. that the circuit outputs 'yes' on a good proof for a 'yes' instance with probability pu = ^■ 
In the language of spins, the circuit U outputs 'yes' when we measure spin up on the output qubit. Therefore, the 
expected value of measuring (Jz on output qubit of the circuit U is bounded from below by 

&vir' > 1 X + (-1) x{l-pu)^2pu-l = \ (34) 
when X £ iyes- Analogously, for x € Lno, it is bounded from above by 

{&:r''<-'^pu + i = -l- (35) 

To solve BQP problems with our automaton, we need to distinguish the 'yes' from the 'no' cases, i.e. we need to 
show that the expectation value of measuring cr*^^^ on the output qubit of our automaton at a random time r < tiq 
is greater than zero in the 'yes' case, and smaller than zero in the 'no' case. The probability of finding a state where 
the computation is finished is piQ. Let us consider a 'yes' instance (x G Ly^s)- Using (PT|) and recalling that the two 
terms in it are orthogonal, we have 

(^S).- = (^^"'(r)| \v'"'Hr)) - (1 -Pio) {^fau{T)\ \vfau{r)) + pio {0u\ <J^l m^^rk ■ (36) 

" V ' 

output of U 

The second term is the circuit output ([34]) . therefore 

{Ou\ <jI:i \0u)^,rk - {'^ii)y:r' > -^pu - 1. (37) 

The first term can be bounded from below (adversarially, i.e. for every time the computation is not finished, the 
output qubit gives the opposite of the correct answer) by 

(^VSl^')>-l- (38) 

Putting it together, we have 

> pwi2pu - 1) - (1 - Pio) = 2pwPu - 1- (39) 



Analogously, for the x ^ L case, we obtain 



-2pioPu + 1. (40) 



Therefore, equations ([391) and ([401) read 



j^^) ) > --0 

'wolves ^ „ >^ 



9 \\ogM J ' 



1 „ / 1 



ir2)„o < + T—TT • (41) 



9 \ log M ^ 

Therefore, we can recognize any language in BQP using the HQCA we described above. 

As an aside, note that there is a way to determine that we obtained a state in which the computation has been 
done with certainty (and thus getting rid of the first term in ([36| ) . We could have chosen to measure all the program 
qudits to the right of the first work qubit and check whether all the S and W are gone. This happens with the above 
probability pio, and the postselected state of the work qubits now surely contains the output of the circuit U . Note 
also that we can think of the state of all the work qubits as the circuit output, as compared to only the last work 
qubit. Nevertheless, thinking only about the last work qubit is enough to ensure universality of our HQCA for the 
class BQP. 

Concurrently with our preprint. Chase and Landahl [23 | found a d = 8 Hamiltonian computer construction uni- 
versal for BQP. However, their Hamiltonian is not translationally invariant, and thus not a HQCA. If we release the 
translational invariance in our c? = 10 HQCA model, we can be sure where the computational qubits are and use the 
gate set {S, SW} instead of {S, W^I} in our program register. This brings the required dimensionality of our model 
to d = 8 as well. However, it remains an open question whether this can be done with keeping the translational 
invariance of the Hamiltonian, i.e. whether a d < 10 HQCA in ID exists. 
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III. THE d = 20 HQCA 



We now present our second construction, a HQCA for a chain of 20-dimensional qudits. As in Section FlI Al we 
describe an encoding of the progression of a quantum circuit U into a set of states of a qudit chain. However, the 
geometry of this set of states \ipt) will be now much simpler, as we can label them by a "time" label t = 1, . . . , T with 
T = poly{N), thinking of the set of states as a "line". The Hamiltonian H20 we construct induces a quantum walk 
on this "line" of states. We conclude by proving that when we let the initial state \ipo) evolve with H20 for a time t 
chosen uniformly at random between and T20 = 0{T\ogT), we can read out the result of the quantum computation 

U with probability P20 > | — O ( j by measuring two of the qudits in the computational basis. 



A. The Construction 



We encode the progression of a quantum circuit U in the form ([2]) (see also Figure [T]) into a set of states jV't) of a 
qudit chain with length L = {2K — 1){N + 1) + 2. As in Section fH Ai each qudit consists of a program register and 
a data register. The data register is again two-dimensional, but the program register can now be in the following 10 
states: 



W, S, I 
► 



the program sequence, 

marked characters in the program sequence, used to propagate 

the active spot to the front (left) of the program sequence, 

apply gate symbol, 

shift program forward, 

a turn-around symbol, 

empty spot (before/after the program). 



Similarly to ^ and (jH), the initial product state 1-00) 
the circuit in Figure [IJi) 



®^=i i\Pj) ^ is given by (we write an example for 



j 


1 
















L 


Pj 




/ 


w 


S 


I 


I 


s 


w 










1 1 wi 


W2 




1 











1 



(42) 



In general, the data register contains N work qubits (labeled w„ in our example) at positions {K — 1){N -I- 1) + 2 + n 
for n = 1 : N (counting from the left). Qubit wat is the designed output qubit for the computation, i.e. once the 
computation is done, wn contains the output of U. Next, the data register contains qubits in the state |1) at positions 
{k — 1){N -I- 1) -I- 2 for k ~ 1 . . . 2K and qubits in the state |0) everywhere else. The I's serve as sequence boundary 
markers. The program register has empty symbols . on the left, and then it contains the program in the form 



/ f/i,i...C/i,jv-i// U2,i...U2,N-iII ■■■II Uk^...Uk,N^1, 



(43) 



gate sequence 



2"*^ gate sequence 



last gate sequence 



with the program written from left to right. In our example (|^^ . the first gate sequence (see Figure [T^) is 11^5' and 
the second gate sequence is SW. Finally, the last qudit in the program register is in the state O , marking an active 
spot in the computation. 

We now give the rules to obtain the sequence of states \ipt) from Ii/jq)- These rules are constructed so that there is 
always only one of them that can be applied to a given state \tpt), thus giving us a unique state |V't+i)- (Also, using 
the rules backwards, one obtains a unique \tpt-i) from lipt))- The first three are 



1 



A 







@ 








A 


® 


— ). 


@ 


B 








® 







A 



(44) 



where A, B stands for either S or /. These rules ensure the passing of the active spot from the back end (right side) 
of the program to the front (left side), without modifying the data register or the order of the gates in the program 
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sequence. Next, we have 



4a 










1 



46 
















[> 








(45) 



After the active spot has moved to the front of the program, there are two possibilities. The turn symbol O can 
change to the apply gate symbol ► (rule 4a), or to the shift program symbol > (rule 4b), depending on whether the 
data qubit below contains the sequence boundary marker state 1. Afterward, for the states containing the apply gate 
symbol we have: 



5a 



► 


A 


X 


y 



A 


► 


A{x,y) 



6a 



► 




1 










1 





(46) 



(Note that rule 6a at the right end of the chain involves only the two particles directly above each other, as no particle 
to the right of them exists.) When applying rule 5a, the apply gate symbol ► moves to the right, while a gate from 
the program sequence is applied to the qubits in the data register below. Applying the rule repeatedly, the ► symbol 
moves to the right end of the program sequence. As an example, we now write out the state IV'12) that we obtained 
from the state \tpo) applying rules 1, 2 (6 times), 3, 4a and 5a (3 times) from the state |'0o)- 



I W S 

1 1 I ...61. 



I I S W . 

1 1 



(47) 



where \. . .9 . . .) stands for the state of the three work qubits after the gates W12 and then were applied to them. 
Let us have a closer look at the marker qubits (all qubits in the data register except for the work qubits w„) and the 
application of rule 5a. The marker qubits stay unchanged for all \^(Jt)■ The gate applied to pairs |0) |1) and |1) |0) of 
marker qubits or the pairs of qubits |1) \qi) and {qn) |1) (the left and right ends of the work qubit sequence) is always 
J, because of the identity gates we inserted between sequences of gates in the program (j43p . Finally, the qubit pairs 
|0) |0) between the 1 markers do not change under the swap operation or the W gate (a controlled gate). 

After the apply gate ► symbol gets to the end of the sequence, it changes into the turn symbol O via rule 6a. 
Note that the boundary markers in the data register are spaced in such a way, that the ► symbol will arrive at the 
right end of the sequence when the qubit below is in the state 1. Using rule 6a, ► will then change into the turn 
symbol O ■ After applying rules 1, 2 (6 times) and 3, the active spot again moves to the left of the program. Because 
the O symbol is now above a marker qubit, rule 4b can be used, and we get a state with the shift program symbol 
O. Finally, here are the last two rules: 



5b 



> 


A 




A 


t> 



6b : 
























(48) 



where again A stands for either W, S or I . Rule 5b makes the program shift to the left while the l> symbol moves to 
the right. Finally, rule 6b deals with what happens when the l> symbol arrives at the end of the program sequence. 
Because of the way we constructed the data register in |'0o), the data qubit below the > symbol will then be in the 
state 0, so that the [> symbol changes to the turn symbol O . The reason why we need to look at the qubit in the 
data register below the O symbol in rules 6a and 6b is that when we apply the rules backwards (making |V't-i) from 
\ipt)), again only one of them applies for each \tpt)- 

After applying rule 1, 2a (6 times), 3 and 4b, the [> symbol appears again and starts shifting the program further 
to the left. After several rounds of this, when the program shifts to the left by iV + 1, rule 4a can be used again 
(as the O symbol will be above a 1 marker qubit), and subsequently, the ► symbol facilitates the application of the 
second sequence of gates to the work qubits. 

After many applications of the above rules, we arrive at the state IV't); for which none of our (forward) rules apply. 



IV't) 



O I W S I I s w 
1 1 I ... 61' . 



1 1 



(49) 



This is the state in which the program has moved to the left of the qudit chain, and all sequences of gates have been 
applied to the qubits in the data register. The state \. . .6' . . .) is thus the output state of the circuit U and the last 
of the work qubits (wn) holds the output of the quantum computation. 

Starting from (gl]), we have constructed the set of states \ipt) for t = . . .T with T = 0{K'^N'^) = poly{N). As t 
grows, these states encode the progress of a quantum circuit U . Let us now think of the geometry of this set of states. 
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They are labeled by a discrete label t, with the state \tpt) obtainable only from the states \ipt-i) and \ipt+i) using the 
above rules and their backward applications. Therefore, the states ji/'t) can be thought of as position basis states on 
a line of length T + 1 

IV't) |t)z»e> (50) 

where t = . . . T. 

Let us choose a Hamiltonian H20 for this system as a sum of translationally invariant terms: 



L-l 6b 

where the terms Pk correspond to the rules l-6b (H ^ . ([¥5 )1 . and ([¥5)1 and act on two neighboring qudits as 



1^—1 DO 



Pi = El® 0|„„®Iaa. (52) 

P2 = ^ (53) 

P3 = E l« ^)( ■@U,p,®Id.,d., (54) 

and 

^4a - I . ►)( . 0|p^,p,®lrf,®|l)(l|,,, (55) 

^^4;, = I . >) ( . 0|p^,p^®lrf, ®|0)(0|,^, (56) 

P5a = E |A ►>(► (57) 
Ag{W,S,7} 

Ps, - 5] |A >)(> ^U,p,®IIrfi,d., (58) 
Ae{w,sj} 

Pea ^ \0 .)(► ■U,p.®|l>(l|d,®Id., (59) 

P6b = \0 .)(> -U.p,® |0)(0|,^®lrf,. (60) 

When thinking of the set of states \ipt) as the set of positions of a particle on a line (|50p . -ff20 becomes 

T-l 

i?H„e = - E + 1| + |i + 1) (<| ) . (61) 

t=Q 

This is the Hamiltonian of a (continuous-time) quantum walk on a line of length T + 1. Therefore, H20 induces a 
quantum walk on the "line" of states \^pt) of the qudit chain of length L. 



B. Required Evolution Time Analysis 



The final step of our model of computation after initializing the qudit chain in the state |V'o) and evolving with H20 
for time r is to read out the output of the computation. As in Section fll Bl we need to ensure that the probability 
of finding the chain of qudits in a state where the computation was performed completely is high. To raise this 
probability, we choose to pad the program (K sequences of gates) with another 5K sequences of identity gates and 
redo the construction in the previous section. The length of the qudit chain thus becomes L = {2{6K) — 1)(A^-|- 1) + 2. 
The states li^tyT/e) (with T modified) now all contain the result of the quantum circuit U in the readout qubit wn, 
as the relevant gates have been applied to the work qubits in those states. Note that as the extra identity gates pass 
by, the state of the work qubits does not change. 

The readout procedure consists of two steps. First, measure the qudit pl-k{n+i) in the program register (the 
qudit with distance from the right end of the chain equal to the length of the original program). Let us call p2o the 
probability to measure . (which would mean the program has moved to the left of the qudit we just measured). 
When this happens, we are assured we have a state in which the computation was done. Second, we measure wn, 
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the last of the work qubits, and read out the result of the computation U. We will now prove that when we choose 
to measure Pl-k{n+i) ^ random time < r < T20 with T20 = poly{N), the probability P20 of obtaining the state 
. is close to |. 

To simplify the notation, let us label the states j-i/'t) as \t). In this basis, the Hamiltonian ([5T|) is the negative of the 
adjacency matrix of a line graph with T + 1 nodes. For the analysis of time evolution with H we refer the reader to 
Appendix [X] We now use the following lemma about a quantum walk on a line (proved in Appendix [X)) : 

Lemma 2. Consider a continuous time quantum walk on a line of length T+l, where the Hamiltonian is the negative 
of the adjacency matrix for the line. Let the system evolve for a time t chosen uniformly at random between and T20 , 
starting in a position basis state \c) . The probability to measure a state \t) with t > T/6 is then p2o > | — O (^T^) ■ 

This implies that when we initialize the qudit chain in the state lipo) (corresponding to the leftmost state on the 
line |c) = |1)) and let it evolve with H for a random time r < T20 with T20 = 0{T\ogT), the probability to find a 
state with t > T/Q is close to |. Therefore, when we measure the program qudit Pl-k{n+i), we will obtain . with 
probability close to |. Finally, when we subsequently measure the work qubit wn, we will obtain the result of the 
quantum circuit U. 

Note that we can also avoid this postselection procedure and simply measure the output qubit. The analysis of the 
outcome would then follow what we did above in Section fll Bl resulting in (PT|) again, with M replaced by T. 
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APPENDIX A: QUANTUM WALK ON A LINE 

Here we analyze the quantum walk on a line and prove two useful lemmas used in Section IIII Bl and Appendix |B] 
Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is the negative of the 
adjacency matrix for the line 

L-l 

H^ = -Y,{\j){j + l\ + \j + l){j\). (Al) 



The eigenvalues of this Hamiltonian are 

A, = —2 cos I 

L + 1 

for j — 1 . . . L, while the corresponding eigenvectors [(/''■-'•') — X]fe=i '^fe'^ 1^) have components 



Aj = -2cos ( ) , (A2) 



.0) - . sin ( ^ 

L+1 \L+1 



-V^^sin(^). (A3) 



Consider the time evolution of a particular basis state |c). The probability of finding the system in a basis state |m) 
at some time r can be found by expanding |c) and |m) in the basis of the eigenvectors (|A3P : 



Pr{m\c) = \{m\e-^^^\c)\' = ^ e-*(^^-^'^-)-</>(^)#*0«*0« . (A4) 

Because the time evolution (according to the Schrodinger equation) is unitary, this probability pr{m\c) does not 
converge. On the other hand, let us define the time average of p7-(m|c) for time < r < T20 as 

Pr^a{m\c) = — / Pr{m\c)dT. (A5) 

T2Q Jo 
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As we will show below in Lemma 1, this average probability distribution does converge to a limiting distribution 
7r(m|c), defined as the T20 00 limit of the average probability distribution (|A5D . All the eigenvalues (IA2p are 
different, so we can express the limiting distribution as 

L 

7r{m\c)= lim P..o(m|c) = |0« , (A6) 

j = l 

which in our case is 

I > 2 + Sm,c + Sm.L+l-c 

= 2iLTT) ■ 



According to the following lemma, the average probability (jASP converges to the limiting distribution TT(m\c). 

Lemma 1. Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is the negative of 
the adjacency matrix for the line. Let the system evolve for time t < T20 chosen uniformly at random, starting in a 
position basis state |c). The average probability distribution pr2oi'\c) converges to the limiting probability distribution 
7r(-|c) as 

J2 \Pr,oHc) - Am\c)\ < O (—) . (A8) 



m — l 



T20 



Proof. To prove our Lemma 1, we recall Lemma 4.3 of [ij for the total variation distance of the probability distribution 
pT2o from the limiting distribution, saying 

^|p.,„(m|c)-7r(m|c)| < — ^ ' ' (A9) 

Using (|A2|) and (|A3p . we can bound the expression on the right of (|A9p . When j is close to k, i.e. \j — A;| < Ci, we 
can obtain 



a(j)|2 



On the other hand, for \ j -~ k\ > Ci we can write 



< 2. (AlO) 



|A, -Afcl L + 1' 
with Ci and C2 constants independent of L. Inserting into (|A9p . we have 



< (All) 



Y: \PrMc) - nim\c)\ < ^ + ^ = O f Ay (A12) 
^1 T20 T20 \T20J 

which concludes the proof. □ 

Using Lemma 1, we will now prove a useful result utilized in the time analysis of the d = 20 HQCA in Section (jllip . 

Lemma 2. Consider a continuous time quantum walk on a line of length L, where the Hamiltonian is the negative 
of the adjacency matrix for the line. Let the system evolve for a time r < T20 chosen uniformly at random, starting 
in a position basis state |c). The probability to measure a state \t) with t > L/6 is then bounded from below as 

Proof. The probability to measure a state |t) with t > L/6 at time t < T20 chosen uniformly at random is 

P20 = Y Pt2o{Mc)- (A13) 
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Starting with (jASp . we have 



O ( — 1 > 

^20 



\Pt2oMc} - n{m\c)\ 

ni—l 



L 



> 



5 



(AM) 
(A15) 

(A16) 
(A17) 



Therefore, the probabiHty of finding the chain in state \tpt>L/6) & random time r < T20 is thus bounded from below 

by 



5 L 

6 \T20 



(A18) 
□ 



APPENDIX B: DIFFUSION OF FERMIONS ON A LINE 



We now prove Lemma 3, a result about the mixing of a discrete free fermion gas. 
Lemma 3. Consider the state 



1*0) 



"/M+2M 



|0) 



(Bl) 



of 2M fermions on the right end of a line with L — {f + 2)M sites. Let the system evolve for a time chosen uniformly 
at random between and tiq with the Hamiltonian 



L-1 

Hf^-Y ^h+i + h.c. 



(B2) 



and measure the number of fermions in the region 1 < x < fM. The probability to measure a number greater than 
Mz.pio>f^-o(4). 

Proof. Let us start with the outline of the proof. We look at the fermionic system in both first and second quantization 
to obtain an expression for the time evolution of the creation and annihilation operators in the Heisenbcrg picture, 
mapping it to a quantum walk on a line. We then consider the observable X, the number of particles sufhciently 
far from the right end of the line. We will show that when we choose the time to measure X uniformly at random 
between and tiq, the expected value we will obtain is approaching a number close to 2M. To show this, we will 
express the expected value of X in the time-averaged state of the system using the results from a quantum walk on 
a line. Finally, because the number of particles in the system is 2M, we will deduce that the probability to measure 
a number less than M is then small. 

Observe that Hf is the Hamiltonian of a free fermion gas on a line in second quantization (a special case of the XY 
model). The time evolution of the state |^'o) can be obtained by looking at the problem back in the first quantization, 
where we write |^'o) as 



1*0) - 



\(l)fM+l) <S) \(f>fM+2) • • • (Xi \4>fM+2M) 



(B3) 



with = \j) in the position basis and [•] the standard antisymmetrization operator. We first solve for the time 
evolution of the corresponding one-particle wavefunction \4>j{T)) with the Hamiltonian 



L-1 



+ 01), 



(B4) 
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and then obtain the solutions for the many-particle problem by antisymmetrization as 

|*(r)) = |0/M+i(r)) «) |(/>/m+2(t)) (g) •• • (g) |(/>/M+2Af(T)) . (B5) 

The eigenfunctions of Hi (quantum walk on a line) are plain waves (as in (jA2p and (jA3p ). and the time evolved states 
|0j(T)) thus readily available. Let us define the unitary matrix m(t) by 

L 

b{r)}=Y,uMr)\k}. (B6) 

fc=i 

Returning to the second quantized system, the time evolution of the creation and annihilation operators in the 
Heisenberg picture is then 



k=l 
L 



k=l 

Consider now the observable X, the number of particles in the first fM sites of the line with length L = {f + 2)M 

fM 

X = ^n™. (B8) 

m—l 

Its expectation value at time r is 

fM 

Er{X) ^ ^(vl/(r)|n„jvl,(r)). (B9) 

m—l 

The number operator for site m is hm = bl^bm. We can go to the Heisenberg picture and use (|B7[1 to write 

(vl'(r)|n™|vl/(r)> = (^ol W l*o> (BIO) 

L L 



J2 E ^™c(TX,(r) {^o\ h% 1*0) (Bll) 

c=l d=l 

L 

Y,\u^,{Tf {^^\b%\^a) (B12) 
1 

L 

E Wrnc{T)\\ (B13) 



c=l 

L 



c=/M+l 

where each term |umc(''")|^ — Pt{iti\c) can be thought of as the probability of finding a particle at site m at time r 



when it started from the site c and performed a quantum walk on a line, according to (|B4p . Inserting this into (jB9[) . 
the expected number of particles not in the rightmost part of the chain at time r is 

L I fM \ 
Er{X)^ E E^'-Hc) . (B14) 

c=fM+l \m=l / 

Let us now choose the time r uniformly at random between and rio- The average value of X (the expectation 
value in the time-average state) is 

EruXX) = — / Er{X) dr. (B15) 



no Jo 

For a quantum walk on a line, the time-averaged probability (|A5[) of finding a particle that started at position c 
at final position m converges to the limiting distribution (|A7p according to Lemma 1 (|A8[) proved in Appendix [X] 
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Using this fact, we can show that the expectation value Erig{X) in the time-averaged state converges to the hmiting 
expectation value 

m<fMc>fM 

as 

\Er,,XX)-EiX)\<o(^^y (B17) 
Recalling the limiting probability distribution for a quantum walk on a line of length L (|A7p , we have 

^(^) = E E (B18) 

m<fM OfM 

= fM X 2Af X — ^ — - + 2M X — - (B19) 
2(L + 1) 2(L + 1) ^ ^ 

= 2M(^y^) +0(1). (B20) 

Putting this into ([BTtI) . the average value of X when the time t < tiq is chosen uniformly at random is bounded 
from below as 

E.JX)>2M{jL.)-o{i^)^ (B21) 

We want to find the probability of measuring X > M . First, the maximum possible value we could measure at any 
time is 2M, the number of particles in the system. Second, the average value Er-^g{X) at time r chosen randomly is 
close to 2M. Therefore, the fraction A of times at which we measure a number significantly lower than 2M must be 
small. Let us bound A in the worst case scenario. This is when each unsuccessful measurement yields X = M , and 
each successful measurement gives us 2M. We then have 

AM + (1 - A)2Af > ^rio(^): (B22) 
A < — Zl^l}_J. . (B23 

Hence we arrive at the desired bound on the probability to measure X > M: 



□ 
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